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Abstract 

We consider a rotating Bose-Einstein condensate in a square optical lattice 
in the regime in which the Hamiltonian of the system can be mapped onto 
a Josephson junction array. In an approximate scheme where the couplings 
are assumed uniform, the ground state energy is formulated in terms of the 
vortex configuration. The results are compared with experimental results 
and also previously reported results for frustrated XY model. We also show 
that vortex configuration is robust with respect to change of couplings and 
therefore the results remain valid when we consider more realistic model with 
non-uniform couplings. 
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1. Introduction 



After the experimental realization of Bose-Einstein condensation (BEC) [l[ , 
it becomes an interesting topic in science ji], 0, 0, Isf. This physical system 
has close relation with other important physical problems in condense mat- 
ter, e.g. Bloch Oscillations Ig, j, [sl, Wannier-Stark Ladders[10l]. Josephson 



junction arrays 



transition 18, 19, 20, 21 
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17[ and superfluid to the Mott-insulator 
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Behavior of a rotating BEC in an optical lattice is similar to that of a 



Superconductor in the magnetic field [23|, l2J, |25|. When there is just trap 



potential without the lattice structure, the Abrikasov vortex lattice can be 



observed just like type II Superconductors [ll|, ll2|, [IJ, ll5|, l26|, |27|, l28|, [29 



30, iMl- When the lattice structure is added to system, the question of 
the vortex lattice in ground state becomes more complicated 32] . There are 
many experimental studies on this system with optical lattices of different 

With some criteria 



symmetries, specially with square lattice [33|, [34, 



35 



this problem can be mapped to problem Josephson Junction arrays (JJAs) 
and numerical experiments show that vortex lattice is same for both the 
systems 



12, 14 



We focus on study of problem in this regime and extend our previous 
results [io*] for an optical lattice with square symmetry. We calculate vortex 
configuration for few rational rotation numbers and show that there is good 
agreement with experimental results on rotating BEC and reported result 
for Josephson junctions lattices |37| . 

Structure of paper is as follows: In section II we formulate problem in 
JJA regime. In section III we focus on vortex lattice with rational rotation 
numbers and compare our result with previous experimental and numerical 
results. In the last section main results are highlighted and suggestions for 
future works have been given. 



2. Rotating BEC in an optical lattice 



Hamiltonian for a BEC in an optical lattice with a rotation frequency Qz 



IS, 



H 



-ihV - MVLz X r) 
2M 



+ + (1) 



where M is the atomic mass and g = Airh^a/M the coupling constant with 
the s-wave scattering length a. Conservation of the total number of particles 
is ensured by chemical potential /i. The external potential Vext = Vh + Vo 
consists of two parts: modified harmonic potential Vh = Mlu^ — n^)r^/2 + 
Mu^z^ /2 and the potential of the optical lattice Vq, which may be chosen as 

For example, the potential 
Vo[sin^ (kx) + sin^{ky)], with 



periodic |38[, quasiperiodic |38l. |39| . or random |40 
with square symmetry can be written as Vq 
the periodicity Tc/k. 

For a large Uz, we can suppose that the system is frozen in the axial 
direction. If the energy due to the interaction and the rotation is small 



2 



compared to the energy separation between the lowest and the first excited 



band, the particles are confined to the lowest Wannier orbitals[lJ|. Therefore 
We can write ip in the Wannier basis as 



N 



"^{r) = ^aiWi{r)expC-^ A{r').dr'), (2) 

where A = Qz x r is the analog of the magnetic vector potential, wjlr) = 
w{r — ri) the normalized Wannier wave-function localized in the i-th well and 
rhi = a\ai the number operator. Substituting [2] in the Hamiltonian [1] leads 
to the Bose-Hubbard model in the rotating frame: 

^ = -t ^ (ajaj-e^^*^ + h.c.) 

<i,j> 

+ ^(ci -/i)mi + ^^mi(mi - 1), (3) 

i i 

where means i and j are nearest neighbors, t — j drw*{—h'^'V'^ /2m + 
Vo)wj is the hopping matrix element, €i = J drw*VhWi an energy offset of 
each lattice site and U = g J dr\wi\'^ the on-site energy. The rotation effect is 
described by Aij proportional to the line integral of A between the i-th and 
j-th sites: Aij = {m/K) f^^ A{r').dr'. 

If the number of atoms per site is large (n^ ^ 1), the operator can be 
expressed in terms of its amplitude and phase, and the amplitude by the c- 
number as a, ~ ^Jrnit^^K Using rhi = nii — id/d9i and 9i = 9i, Eq. [3] reduces 
to the quantum phase model: 



<ij> 



where Jij = 2t^mimj. Also, the atoms are assumed to be distributed such 
as they satisfy the condition ej + Uni = fi. The magnitude of Jij decreases 
from the central sites outward since rrii has the profile of an inverted parabola 
as we can see from solution of ID case in Fig. [H 

Eq. 2.4 is just the Hamiltonian of a lattice of small Josephson junctions 
with inhomogeneous coupling constants Jij. With Jij ^ U, we can neglect 
the kinetic energy term; Moreover, we suppose that coupling constants are 
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Figure 1: Solution of GP equation for a periodic one dimensional potential. It can be seen 
that l^'p is enveloped by an inverted parabola shape. 



equal Jij = J. We also check the robustness of the results for non-uniform 
coupling. 

This regime can be achieved when the depth of the optical lattice Vq is 
between 18 and 25Eji [E^ = h'^k'^/2M is the recoil energy) and the average 
particle number ITojill. Then we have: 



= - ^ cos(7i 



(5) 



<i,j> 



with = 6i — 6j — Aij. Our purpose is to find the minimum of this Hamilto- 
nian as a function of the condensate phases 6i with two constrains as follows. 
First constraint originates in the single valuedness of the wave function which 
leads to the quantization of vorticity: 



y^7»i = 27i{nk - akf) 



(6) 



where sum is over A;th plaquette, is the area of kth plaquette, and is in- 
teger. In the case of a condensate with the rigid-body rotation, the frustration 
parameter is defined as / = 2Q/k with the quantum circulation k = h/m. 
In two dimensional arrays, frustration represents density of vortices per the 
number of the plaquettes j37l |. 

The second constraint comes from conservation of the particles numbers. 
The particle current from the nodes can be assumed as the derivation of the 
energy function with respect to the corresponding phase 421] . This leads to 



0, 



(7) 
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where sum is over all j which is connected to ith node. 

Since Hamiltonian is periodic, can be [a^/J or [a^/J + 1, and we can 
replace them by Uk which is zero or one, and replace akf with akf — [flfc/J- 
These two set of equations can completely determine all jij for a given set of 
Uk- Also, for a lattice, number ofplaquettes plus number of nodes is equal to 
number of connections minus onqj; therefore one of the equations (e.g. one of 
particle conservation equations) can be neglected. Assuming the particle flow 
from the nodes are small, we can linearize the EqJ7]i.e. sin'jij fa 7jj. Now, we 
have a linear set of equations which must be solved for a given set of integers 
Uk, i.e. for a vortex configuration. Energy of this vortex configuration is 
denoted by E^n}- Then instead of original minimization scheme, we can 
use the set of these energies for minimization of Hamiltonian. It means 
that minimum of Hamiltonian correspond to a vortex configuration which 
minimizes these set of energies, with number of vortex configuration equal 
to 2-'^ where N is number of plaquettes in the lattice. 



3. vortex configuration in ground state: ciieckerboard structure 



We consider the square lattice and present the result of above formulation 
for some rational values of /. When the number of the plaquettes becomes 
large, the checkerboard vortex structure of ground state can be observed |3 71 . 
The same pattern can be observed in the trapped BEC[33|, [SJ, l35|. We show 
that our formulation has the same result. 

Set of the equations Eland the linearized version of[7]for the square lattice 
can be written as 



C7 = b, (8) 

where 7 is a the vector representation of the gauge invariant phase differences 
and C is the matrix of coefficients:, 

Ci,k = -1 (9) 

Q,k-i = 1 (10) 

Ci^k-m = 1 (11) 

Ci,fc-m-l = — Ij (12) 



This is Euler formula for polygon lattice excluding exterior face. 
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for i = 1, ..,nm — 1 where k = \i/m\ + i, n is number of rows and m is 
number of columns; and for i = nm, 2mn — m — n, 





= 1 




= 1 


Ci^k+m—l 


= -1 


Ci^k—m 


= 1, 



(13) 
(14) 
(15) 
(16) 

where k = {m — 1) [{i — mn + l)/mj . 6 is zero for i = 1, ...,nm — 1 and 

bi = 2TT{n{i — mn + 1) — a{i — mn + 1))) (17) 

for i = mn, 2mn — m — n. These equations can solved for any given m, n 
and set of and Oj. The result for the minimum vortex configuration for 
an square lattice with = 1, m = = 4 is shown in Fig. [2] for / = 1/2. 
In this case for < / < 1/2, rii could be zero or one. The plaquettes with 
points in Fig. [2] are those with rij = 1. The minimum vortex configuration 



shows the checkerboard pattern which is known for frustrated XY mo del [37 



For / = 1/3, the minimum vortex configuration, is shown in Figj3l Again 




Figure 2: Minimum vortex configuration for 4x4 square lattice for / = 1/2 
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the checkerboard pattern is seen for the vortex configuration in the ground 
state. In both cases, we can see that there is an extra vortex in the boundary 
of the lattice. This is a boundary effect and can be ignored in large lattice 
limit. These results are in agreement with experimental results on rotating 



BECl33|, M, l35 



Figure 3: Minimum vortex configuration for 4x4 square lattice for / = 1/3 



When number of plaquettes increases, number of possible vortex con- 
figurations, increases and we can use a combination of this method with 
Monte-Carlo method 36 . 



4. Robustness of minimum vortex configuration with respect to 
coupling constants 

We see that in our approximation, the coupling of JJA changes as the 
inverted parabola. Here, we discuss the effect of changing the coupling 
constants on the minimum vortex configuration. We consider an inverted 
parabola mj, 

m,j = 16u(^ ~Nr- -Nc- 1)/[{Nc + l)\Nn + 1)^] (18) 
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for ijth site of square lattice, where Nc is the number of columns and 
is the number of rows. The coupling between two nearest neighbor sites 
on the lattice can be found from J{ij)^{k,i) = y/^ij^ki which gives a better 
approximation comparing to the uniform couplings. Result for the same 
situation as Fig. |2l is shown in Fig. |H It is seen that the minimum vortex 
configuration is not affected by the change of the coupling configuration. 
We have also checked the results for random mij and we found the same 
minimum vortex configuration. 




Figure 4: Minimum vortex configuration for 4 x 4 square lattice for / = 1/2 with non- 
uniform coupling. 

5. Conclusion 

We study a rotating BEC in a square optical lattice in a regime which 
Hamiltonian of the system can be mapped onto JJA. In this regime, we 
formulate ground state of system in terms vortex configuration. Our results 
in uniform coupling case, show that vortex configuration in ground state has 
checkerboard pattern. We have also checked the non-uniform coupling in the 
form DJD, which D is a diagonal matrix, and J is original uniform coupling 
matrix. The results for this case show that minimum vortex configuration is 
robust with respect to the change of coupling matrix. 
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